This paper presents a numerical method for a class of singularly perturbed parabolic partial differential equations with integral boundary conditions (IBC). The solution to the considered problem exhibits pronounced boundary layers on both the left and right sides of the spatial domain. To address this challenging problem, we propose the use of the implicit Euler method for time discretization and a finite difference method on a well-designed piecewise uniform Shishkin mesh for spatial discretization. The integral boundary condition is approximated using Simpson's $\frac{1}{3}$ rule. The presented method demonstrates almost second-order uniform convergence in the discretization of the spatial derivative and first-order convergence in the discretization of the time derivative. To validate the applicability and accuracy of the proposed method, two illustrative examples are employed. The computational results not only accurately reflect the theoretical estimations but also highlight the method's effectiveness in capturing the intricate features of singularly perturbed parabolic partial differential equations with integral boundary conditions.

In this paper, an approximate solution of a nonlinear parabolic partial differential equation is obtained for a non-uniform mesh. The scheme for partial differential equation subject to Neumann boundary conditions is based on cubic B-spline collocation method. Modi, ed cubic B-splines are proposed over non-uniform mesh to deal with the Dirichlet boundary conditions. This scheme produces a system of , rst order ordinary differential equations. This system is solved by Crank Nicholson method. The stability is also discussed using Von Neumann stability analysis. The accuracy and efficiency of the scheme are shown by numerical experiments. We have compared the approximate solutions with that in the literature.

In this paper, a reliable numerical scheme is developed and reviewed in order to obtain approximate solution of time fractional parabolic partial differential equations. The introduced scheme is based on Legendre tau spectral approximation and the time fractional derivative is employed in the Caputo sense. The L2 convergence analysis of the numerical method is analysed. Numerical results for different examples are examined to verify the accuracy of spectral method and justi, cation the theoretical analysis, and to compare with other existing methods in the literatures.

We deal with the initial-boundary value problem for a quasilinear degenerate parabolic equation with inhomogeneous density and absorption, which appears in a number of applications to describe the evolution of diffusion processes, in particular non-Newtonian flow in a porous medium. We discuss the extinction of solution and the finite speed of propagation of perturbations.

In this paper, we solve the Caputo's fractional parabolic par-tial integro-di, erential equations (FPPI-DEs) by Gaussian-radial basis functions (G-RBFs) method. The main idea for solving these equations is based on the radial basis functions (RBFs) which also provides ap-proaches to higher dimensional spaces. In the suggested method, FPPI-DEs are reduced to nonlinear algebraic systems. We propose to apply the collocation scheme using G-RBFs to approximate the solutions of FPPI-DEs. Numerical examples are provided to show the convenience of the numerical scheme based on the G-RBFs. The results reveal that the presented method is very efficient and convenient for solving such problems.

In this work, the inverse quasi-linear pseudo-parabolic problem was investigated. We demonstrated the solution by the Fourier approximation. The inverse problem was first examined by linearizing and then used implicit finite difference schema for the numerical solution.

In this paper, a variational iteration method (VIM), which is a well-known method for solving nonlinear equations, has been employed to solve an inverse parabolic partial differential equation. Inverse problems in partial differential equations can be used to model many real problems in engineering and other physical sciences. The VIM is to construct correction functional using general Lagrange multipliers identified optimally via the variational theory. This method provides a sequence of function which converges to the exact solution of the problem. This technique does not require any discretization, linearization or small perturbations and therefore reduces the numerical computations a lot. Numerical examples are examined to show the efficiency of the technique.